Ultrapure Multilayer Graphene in Bromine Intercalated Graphite 
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We investigate the optical properties of bromine intercalated highly orientated pyrolytic graphite 
(Br-HOPG) and provide a novel interpretation of the data. We observe new absorption features 
below 620 meV which are absent in the absorption spectrum of graphite. Comparing our results 
with those of theoretical studies on graphite, single and bilayer graphene as well as recent optical 
studies of multilayer graphene, we conclude that Br-HOPG contains the signatures of ultrapure 
bilayer, single layer graphene, and graphite. The observed supermetallic conductivity of Br-HOPG 
is identified with the presence of very high mobility 121,000 cm^V~^s~^ at room temperature 
and at very high carrier density) multilayer graphene components in our sample. This could provide 
a new avenue for single and multilayer graphene research. 



PACS numbers: 72.80.Vp, 73.22.Pr, 78.67.Wj 

Recently several allotropes of carbon including 
Buckminster-fullerenes, carbon nanotubes, and graphene 
have been discovered and have attracted much 
attention[lHZ]- These carbon-based materials have 
unique crystal structures and interesting and useful phys- 
ical, chemical, and mechanical properties, leading to po- 
tential applications in various fields including carbon 
based electronics. Multilayer graphene has become es- 
pecially important because of the known unique Dirac 
cone band structure in single layer graphene with its chi- 
rality, no n- trivial Berry phase, half integral Hall plateaus 
etc pHTJ. The parent material for graphene is graphite, 
which has a layered structure. Each layer of graphite 
consists only of carbon with hexagonal lattice structure. 
The in-plane carbon-carbon interaction provides a strong 
covalent bond and inter-plane interaction is a weak van 
der Waals force. One can easily insert atoms or molecules 
between the sheets in graphite through intercalation pro- 
cesses. In most cases Li, Ca, and K have been used for 
intercalation jSHTT] . Studies of Br intercalated graphite 
mostly focused on the structural properties of the inter- 
calant layer itself[T2l|T3]. For our considerations here it is 
important to note that the c-axis conductivity at high Br 
concentrations is reduced by more than factor of 5 over 
its value in graphite [H [14] while for K it is increased by 
a factor of ~ 25. Thus Br2 leads to more electrically iso- 
lated graphite planes (two dimensional graphene) while K 
integrates them into a more three dimensional electronic 
structure. 

In this study, we intercalated bromine molecules be- 
tween graphite sheets by exposing highly oriented py- 
rolytic graphite (HOPG) to a bromine gas environment. 
We prepared three bromine intercalated samples with dif- 
ferent exposed time, 15 minutes, 30 minutes, and 100 
minutes, respectively, labeled as 15-min Br-HOPG, 30- 
min Br-HOPG, and 100-min Br-HOPG, respectively. In 
general, the intercalation can introduce charge carriers in 
the graphite sheets and can make the distance between 



sheets larger. However, it is not clear how uniformly Br 
atoms are distributed between the graphite layers and 
how much charge transfer occurs. It is known that Br- 
HOPG shows carrier density increases, supermetallic be- 
havior, and reduced scattering rate at high intercalation 
levels [14l[T5]. The reduction of the scattering rate at high 
Br intercalation is yet to be understood. Here, we focus 
on this issue and on additional new features, which ap- 
pear below 620 meV in the optical spectra of Br-HOPG. 

In Fig. 1 (a) we show the in-plane optical con- 
ductivities of four samples including bare HOPG, ob- 
tained from the measured reflectance via Kramers-Kronig 
analysis [16]. The optical conductivity of the HOPG 
at room temperature agrees well with that previously 
reported[17, 18 . The optical conductivities of the Br- 
HOPG samples show dramatic changes as compared with 
HOPG. Some of the properties of Br-HOPG have al- 
ready been reported, such as the supermetallicity, aris- 
ing from a high Drude plasma frequency and a scattering 
rate reduction [TT. The spectra show the appearance of 
new absorption features in mid-infrared region (marked 
with arrows in the figure) and an overall suppression of 
the conductivity in the spectral region from mid-infrared 
through near infrared with increasing Br intercalation. 

To understand the conductivity suppression of Br- 
HOPG in mid-infrared and near infrared region we take 
the optical conductivity of bare HOPG near 620 meV 
(the relatively flat part of the spectrum) as a reference. 
This reference is the parallel contribution to the conduc- 
tivity of each graphite layer separately [18 . We normalize 
the conductivity at 620 meV to Ggo, the effective opti- 
cal conductance of a single graphite sheet (graphene) in 
HOPG. The normalized conductance can be written as 
G{uj)/Gso = dc cri(cj), where dc is the c-axis lattice con- 
stant and ai{uj) is the real part of the optical conductiv- 
ity. In bare HOPG dc = 3.4 A. In the Br-treated samples, 
the normalized conductance yields the c-axis spacing, 
because it is related to the number density of graphite 
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FIG. 1: (Color online) (a) The room temperature conductiv- 
ity spectra of bare HOPG and three Br-HOPG with different 
intercalation times. We indicate the five differences between 
bare HOPG and Br-HOPG with vertical arrows (see in the 
text), (b) We show a schematic diagram of Br2 intercalated 
graphite: (upper) a top (in-pane) view and (lower) a side view 
(along c-axis) of A-A line cross section. This diagram is based 
on results by Erbil et a/.|13j. 



sheets along the stacking axis. In the inset of Fig. 2 (a) 
we display the extracted c-axis constant, dc as a function 
of (Br-intercalation time)^/^ for all four samples; values 
are scaled by the c-axis lattice constant of bare HOPG. 
The c-axis lattice constant increases as the amount of Br 
intercalants increases, because the Br atoms go between 
the graphite layers and push the graphite sheets further 
apart as sketched in Fig. 1 (b). In the early stage of in- 
tercalation (up to 30 minutes), dc increases as the square 
root of time, as found in Ref. [M] for Br uptake, and then 
begins to saturate reaching dc — 5.2 A for the 100-min 
sample. This spacing is smaller than the 7.0 A Br-layer 
spacing found in structural studies [13] but bigger than 
the average of 4.4 A. Our effective c-axis lattice constant 
is an average spacing. We can not tell whether the Br 
atoms go between every graphite sheet or not. 

It is interesting to compare further the conductance 
of our Br intercalated samples with the conductance of 
HOPG shown in Fig. 1 (a). We emphasize five differ- 
ences: a broader Drude mode, a peak near 85 meV, a 
reduction of spectral weight near 300 meV, a new peak 
near 400 meV, and a valley near 450 meV. The peak 
around 85 meV in the 15-min sample (red dashed line) 
shifts to lower energies and appears to vanish in the 100- 



min sample (blue solid line). The reduction at higher 
photon energies below roughly 300 meV is accompanied 
by increased optical spectral weight in the Drude spec- 
trum. The new peak near u = 400 meV is common to all 
three doped samples. This peak is followed by the valley 
before the conductance rises to its constant uniform back- 
ground at high energies. To gain insight into the possible 
origin of these five features it is useful to recall recent the- 
oretical and experimental studies of the conductivity of 
single [7| [T8H2Q] and double layered graphene [2TH25] . 

For a single layer in the clean limit a remarkably sim- 
ple formula of the conductance, Gs(a;), was obtained 

[la Uni namely %^ = 4|/i| S{uj) + O^f - 

where Gso = 7Te^/{2h) is the conductivity of a single 
graphene sheet, e is the charge on the electron, and h 
is the Plank's constant. When broadening is included, 
the Dirac delta function in the first term is replaced by a 
Drude spectrum. The second term is a Heaviside function 
giving a constant universal background for photon energy 
uj > 2|/i| with /i, the chemical potential describing charg- 
ing. The optical spectral weight under the Drude spec- 
trum which comes from the intraband transitions is ex- 
actly equal to the missing weight under the background, 
which starts at 2|/i| and has its origin in the interband 
transitions from lower to upper Dirac cones. Identifica- 
tion of this onset allows a measure of and hence of 
the charge transfer from the Br to the graphite planes. 

For the bilayer a slightly more complicated formula for 
the conductance, Gb{uj)^ can be derived [25 and is 
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2Gso = e^/2fi, which is twice the conductance of a single 
graphene sheet. The new feature is that now the hop- 
ping parameter between planes, 7, enters. The absorp- 
tion peak (a Delta function in the clean limit) at = 7 is 
a result of the energy splitting between the bands. One 
set of bands is displaced with respect to the other by 
7 and at finite fi this leads to a degeneracy of optical 
transitions between the split bands. For graphite [27H29] 
7 = 400 meV and we identify this new absorption with 
a peak in our data for (Ji{uj) seen in all three Br interca- 
lated samples around 400 meV. As the Br is intercalated 
the system acquires some of the characteristics associated 
with bilayer graphene. 

In Fig. 3 (a) we show numerical results, [Gs{oo)/Gsq + 
Gi){uj) /Gi)q\/2^ based on equations (1) and (2) for an 
equally weighted sum of a single and a doubled layer of 
graphene for several values of the chemical potential. In 
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FIG. 2: (Color online) (a) Room temperature optical conduc- 
tivity normalized at 620 meV, as described in Ref. [18 . This 
normalization allow the extraction of an effective c-axis lattice 
constant, shown in the inset (see in the text), (b) Simulation 
of our data using our model calculation. (See the text for a de- 
tailed description.) In the inset we show a schematic diagram 
of a band structure of bilayer graphene with an anisotropy 
gap, A = 0.447 (not to scale). Here, the gap is caused by a 
built-in potential difference between upper and lower layers. 
The red arrow shows the intraband Drude conductivity, the 
blue the interband transitions, the orange the 7 transition, 
and the dark green the onset of the interband transition. The 
gray dotted arrow shows the transition between the minimum 
gap, Egi. We note that this transition would not appear at 
low temperature or in a highly doped system. 

all cases there is a Drude spectrum at small cj, a con- 
stant universal conductance background at large u and 
a sharp rise in absorption at uj = 2|/i|. There is also a 
peak at = 7, the energy shift between the two sets of 
bands in the bilayer. When twice the chemical potential 
2\fi\ is greater than 7 there is also a region of depressed 
conductivity between 7 and 2|/i|. All these features are 
seen in the blue solid curve of Fig. 2 (a) for the most Br 
doped case and we identify the chemical potential as ~ 
260 meV = 0.667. 

A qualitative understanding of the data for two other 
samples with less Br intercalated requires further elab- 
oration of the results for bilayer graphene. When the 
two layers in the bilayer are not equivalent as would 
be expected if the Br intercalation process is not the 
same on the two surfaces of the bilayer, there can be 
an anisotropy gap, A. Details of the theory can be found 
in reference [25 but here it is sufficient that we model 
the main results. When a gap develops between bands 
that would otherwise be degenerate, part of the Drude 
spectrum will become gapped. We use this idea in the 



results of our own simulations shown in Fig. 2 (b). We 
add the conductance of three subcomponents: graphite 
taken as measured in our experiment and an equal mix- 
ture of single and bilayer graphene. G total {^) / Gtotaio = 
fgGg{uj)/GgoH'i--fg)[Gs{oj)/Gso+Gb{uj)/Gbo]/2, where 
Gtotalo is Gso X the number of total layers, fg is the 
graphite fraction, and Gg/Ggo is the normalized conduc- 
tance of the bare HOPG. We used fg = 0.70, 0.50, and 
0.38 for Br 15-min, 30-min, and 100-min samples, respec- 
tively. These graphite fractions are consistent with the 
values estimated independently using the effective c-axis 
constant 4 = fg 4o + (1 - fg) dcsat (see Fig. 2 (a)), 
where dco is the c-axis constant (3.4 A) of graphite and 
dc-sat is the graphite-bromine-graphite sandwich width 
{c^ 7.0 A)[13 : fg =0.75, 0.63, and 0.48, respectively. 
We also shifted the chemical potential for first and sec- 
ond curves to account for an overall reduction in Br, but 
kept the same ratio between ionized and non-ionized as 
found for the most Br rich sample so the /i's were set to 
O.4I7, 0.487, and 0.657, respectively. Finally for the bi- 
layer component in each case we included an anisotropy 
gap of A = 0.447. To do this, in equation (2) we re- 
place the Drude term A{/j.)d{uj) with a two component 
term fgapA{/j.)S{uj - Egi) + (1 - fgap)A{fi)6{uj) and fgap 
= 0.8, 0.5, and 0.0, for Br 15-min, 30-min, and 100-min 
samples, respectively, where fgap in the portion of gaped 
states and Egi{= ^/S./ (2^ /S? +7^)) [25] is the minimal 
gap between upper and lower bands (see Fig. 2 (b)). 
These percentages roughly represent the idea that as dop- 
ing proceeds the mechanism for anisotropy gap formation 
because less effective as the top and bottom surface of a 
given bilayer should experience more and more the pres- 
ence of the same Br coating. 

The information on the chemical potentials allows us to 
estimate the charge transfer to the graphite sheets. Via 
weight uptake measurements, we know the amount of Br 
intercalated versus the time of exposure to the Br atmo- 
sphere: 2.3, 3.3, and 6.0 at % Br for 15-min, 30-min, 
and 100-min samples, respectively [M]. The bandwidth 
is W ^ 7.0 eV. Because we have a linear band struc- 
ture (Dirac cone) near the Fermi Surface, the number of 
filled states is roughly proportional to where ji is 

the chemical potential. For our 100-min sample, if all 
intercalated Br are ionized the chemical potential would 
be \iJiassumed\ = VO.O6 W = 1.72 cV. The extracted real 
chemical potential from the fit is |/iread = 0.657 — 0-26 
eV, which comes from the ionization. The ionization 
rate is (l/iread/lMassnmedl)^ = 0.023 or 2.3 %. We also 
calculate the ionization rates for other two lower doped 
samples (15-min and 30-min) and get 2.4 % and 2.2 %, 
respectively. The ionization rate can be determined in 
a different way by using the measured change in plasma 
energy [14 and the effective c-axis lattice constant, dc. 
The density of bare HOPG is roughly = 2.16 g/cm^. 
By intercalation the density of carbon will be reduced to 
pm X dcid/dc. For our 100-min sample, the charge carrier 
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FIG. 3: (Color online) (a) Numerical calculations based on 
Eq. (1) and (2) for an equally weighted sum of single and 
bilayer graphene for several chemical potentials. We include 
the temperature broadening in the calculation. The orange 
horizontal arrow shows the degenerate 7-transition. The 
dark green arrow the onset of the interband transition for 
jjL — — I.O7. (b) We show schematic diagrams of band struc- 
tures of both single and bilayer graphene. The red arrow 
shows the intraband Drude transition, the blue the interband 
transitions, the dark green the onset of these transitions, and 
the orange the degenerate 7-transition. 



density is N^r- — 10^^ holes/cm^ p31 and the effective 
c-axis constant is dc = 5.2 A. The density of carbon is 
2.16 g/cm^ X 3.4 A/5.2 A = 1.41 g/cm^, and the corre- 
sponding number of carbon atoms per cm^ is 1.41 g/12.0 
g X Na = 7.08x10^^ cm~^, where Na is the Avogadro's 
number. The total number of intercalated Br atoms is 
Nbt = 0.06 X (7.08 X 10^2) = 4.25 x 10^^ cm-^. The 
ionization rate is Nbt- /^Br — 0.024 or 2.4 %. For other 
two samples (15-min and 30-min) the estimated ioniza- 
tion rates are 3.0 % and 2.6 %, respectively. The two 
approaches give very similar results, and agree with ear- 
lier results. [30] 

We also estimate the conductance of the single layer 
graphene in our samples using the fitting parameters, 
the measured dc resistivities, and the measured carrier 
density[14 : Gsam^ie = fgGg + (1 - fg)[Gs + Gb]/3 ^ 



fgGg + (1 - fg)Gs. So Gs 
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fg). For 100- min sample at room temperature since 
Gsampie = 6.3 X 10"^ [Q-^], Gg = 0.1 X 10"^ [n-^l and 
fg = 0.38, the conductance of the single layer graphene, 
Gs = 1.01 X 10~^ [^~^]- The corresponding mobility, 
Gs/{Ti2d e), is 121,000 cm^V~"'^s~"'^, where n2d is the two 



dimensional carrier density, 5.2x10^^ cm~^ for the 100- 
min sample. The estimated mobility is about 10 times 
larger than that of single layer graphene fabricated on 
a substrate |31j and is comparable to or slightly better 
than that of suspended single layer graphene [32^ even 
though it is at higher temperature (300 K) and at much 
higher carrier density. Interestingly, the graphene in our 
100-min Br-HOPG sample shows higher carrier mobility 
than any reported values; it is " ultrapure!^ graphene. 

In conclusion we investigated bromine-intercalated 
highly-oriented pyrolytic graphite (Br-HOPG) using in- 
frared and optical spectroscopy. New features in the opti- 
cal spectra of Br-HOPG, absent in the optical spectrum 
of HOPG, are observed. The characteristic energies of 
these new features, when compared with published op- 
tical data for graphene, graphite, and especially bilayer 
graphene, allow us to identify them as characteristic of 
the bilayer. Additionally, the chemical potential obtained 
from the simulation tells us that only a few percent of the 
intercalated Br molecules are ionized. The mobility in 
the isolated multilayer graphene components is found be 
very high and provides an explanation for the observed 
supermetallicity of our samples. This is a particularly 
attractive possibility considering the present excitement 
in the field of graphene and related materials [3]-[7]. Mea- 
surements of the optical conductivity of highly interca- 
lated samples under a magnetic field perpendicular to the 
layers would be particularly interesting. We believe that 
our results may open a new avenue to the study of the 
properties of multilayered graphene. 
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